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^ ■ The cosmic ray background consists of several different particles such as pions, 

electrons, protons, etc. They are coming either from space or are generated in 
the high atmosphere. The most numerous generated underground are muons 
that are created when protons and secondary particles interact with atoms in 
the atmosphere. Experiments can be shielded from this cosmic activity by go- 
ing deep underground. In this paper we describe a measurement of the cosmic 
ray muon flux at the Waste Isolation Pilot Plant (WIPP) with an overburden 
of 655 m of rock and salt and located 30 miles east of Carlsbad, New Mex- 



ico (|DOE01f ). Such measurements are useful in assessing the appropriateness 
of WIPP as a site for underground experiments that require shielding from 
cosmic radiation. The depth requirements and adequacy of a particular site 
vary depending on the nature and goals of a particular experiment. This pa- 
per does not attempt to address the appropriateness of the WIPP site as a 
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general purpose undcrgroTind laboratory, an exercise of intense activity within 
the scientific community at present. Such a task is beyond the scope of this pa- 
per. The paper focuses on the measurements of the muon flux at WIPP, which 
serves as an important input in evaluating the impact on specific experiments. 



2 Experimental Setup 

The experiment to measure the muon fiux in the underground of WIPP was 
conducted 655 meters below the surface in the experimental Q-Area. The de- 
tector consisted of two plastic scintillator panels each 305 cm long, 76.2 cm 
wide and 2.54 cm thick. Each panel had two light-guides mounted on the two 
narrow ends. A 12.7 cm diameter photo-tube with an operating voltage of 800 
volts was mounted on each light-guide enabhng a double ended readout. The 
muon detection efficiency as a function of position for each panel was measured 
by placing two small scintillators above and below. If the small detectors trig- 
gered coincident, a muon must have passed through the large panel. The num- 
ber of events in the panel over the total number of triggers gives therefore the 
muon detection efficiency. Each position was tested with 10,000 coincidence 
events. High energy, minimum ionizing muons would deposit about 5 MeV 
when traversing a scintillator paddle at normal incidence. The actual energy 
spectrum deposited in the Patel exhibits a broad tail around this pronounced 
peak due to muons traversing at angles other than normal, straggling effects, 
and the finite energy resolution of the scintillator. Nonetheless, by summing 
the signal heights from the PMT's of each panel with a threshold of 3 MeV, a 
clean muon signal can be obtained well above the electronic noise and low en- 
ergy gamma background associated with a single PMT channel. The detector 
efficiency for muons was found to be 100±1% for all positions. 

For the actual experiment panel 1 was mounted on top of panel 2 in a distance 
of 30.5 cm (see Fig. 1). A fourfold coincidence from the photo-tubes was 
required to record an event. The high voltage on each photo-multiplier tube 
(PMT) was adjusted by demanding an equal pulse height for a ^°Sr calibration 
source place in the center of each panel. The decay product of ^°Sr is 
which decays emitting high energy electrons, which were detected in the panel. 
Fig. 2 shows the electronic diagram. The pre-amplifier signal had a rise time of 
10 ns to 20 ns. It was split up and one signal was fed through a spectroscopic 
amplifier with a shaping time of 0.25 /xs and connected to an LeCroy AD811 
ADC. The other signal was put into a discriminator with a 175 mV threshold 
which was equivalent to an energy of 1.1 MeV. The steep NIM pulses from the 
discriminators were used as input for a coincidence unit. Fig. 3 illustrates the 
timing for the different components during an event. The coincidence window 
had a width of 100 ns. The output of the coincidence was then processed by a 
gate/delay generator and used as strobe for the ADC. The strobe window was 
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chosen to be 0.5 /is wide. This ensured that the peak from the spectroscopic 
amphfier was inside the ADC strobe window. 



3 Data Analysis 



The experiment was conducted during a run-time of 6.17 days. Panel one had 
a slightly better resolution than panel two. This was due to an increased noise 
in photo tube 3. The total charge in a panel was calculated as the sum of the 
outputs of its two PMTs. 

To determine the energy calibration of the detectors Monte-Carlo simulations 
were carried out with Geant 4. The initial energy distribution of the muons 
simulated was derived from the followin g empir ical formula corrected for the 
energy losses in 655 m of rock and salt (|MIY73t ): 



N{P)dP = AP-"" dP with (1) 
a = 0.5483 + 0.3977 In P and 

A = 3.09 X 10"^cm~2s^^sr-^(GeV/c)^\ 

P bei ng in Ge V/c. The intensity distribution of the polar angle 6 was taken 
from (IMIY73I). A simplified angular distribution for medium depth can be 



written as (|MIY73t ) 



9) = 0) cos^-^=^ ^g-8.oxio-4h{scce-i)_ (2) 

with 6) as the flux for a specific angle in cm^^s^^sr^^. Figure 4 shows this 
spectrum at sea level. The energy calibration was done by fitting the maximum 
intensity of the muon peak to the maximum intensity of the Monte-Carlo for 
the detector. A comparison of measured and simulated data is shown in Fig. 5 
for both panels. 

The total muon spectrum for each panel was achieved by adding the ADC 
values of both photo tubes up on an event-by-event base. The different energy 
resolutions of the PMT's were taken into account during the Geant 4 sim- 
ulations (s ee differ ence between upper and lower panel in Fig. 5). For more 
details see fjESCOlh . The shape of the measured spectrum could be reproduced 



within statistical uncertainties. The muon peak with its maximum intensity 
around 4.8 MeV is clearly visible. 

Fig. 6 shows a two dimensional plot of an event by event energy distribution 
for the two panels. The x-axis shows the energy deposition in the lower panel 
the y-axis the energy deposition in the upper panel. Among the events marked 
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with a cross, the muon peak is clearly visible at approximately 5 MeV energy 
deposition in each panel. A cut above 2 MeV energy deposition for each panel 
was chosen to determine the number of muons. The black line represents the 
cutting border. Fig. 7 shows the same cut applied to the Monte-Carlo simu- 
lations. According to the Monte-Carlo simulations 1.24 ± 0.01% of the muon 
events are cut out. The low energy events in Fig 6 are due to PMT noise, 
which are therefore absent in Fig. 7. Table 1 contains the integral information 
of the different regions in the spectrum shown in Fig 6. The errors quoted are 
purely statistical errors. The final muon peak cut is shown by cut 6 and the 
muon number derived from the cut is 



= 5202 ±72, (3) 
where the uncertainty is statistical only. 

The next step is to determine how many counts of the background contributed 
into the signal and how many counts from the signal were lost into the low 
energy background due to the cut. The trigger rate coming from two photo 
tubes in one panel in coincidence was 220±10 Hz. As described above, during 
the experiment the data were taken demanding a quadruple coincidence, which 
basically means a coincidence between the two panels. The random coincidence 
rate can therefore be calculated with 



i/r = 2l/iI/2T, (4) 

where Vr is the random coincidence frequency, Vi and V2 the two frequencies 
and T the coincidence window. With a coincidence window of 100 ± 5 ns the 
rate can be calculated iovr = 9.68 ± 0.79 x 10~^ Hz or 5158 ±43 counts. 
The main part of the background can be attributed to random coincidences 
(see first cut of Table 1) the tail of the background was assumed to be 
exponential and extrapolated into the region of interest. The results are listed 
in Table 2. The number of muon induced events outside the region of interest 
was estimated using the result of the Monte-Carlo simulations (see Fig. 7 and 
Table 1). 

The muon flux can be calculated as follows. 



^total = — , (5) 

where ^totai is the total flux of muons through the panel, Nq the measured 
number of muons, e the efficiency of the detector which includes both the 
physical efficiency of the scintillator and the geometrical efficiency due to the 
setup of the detector and the angular distribution of the muons, a the area 
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covered by the detector (2.33 m^) and trt the run-time of the experiment which 
was 532800 seconds. 



The deviation of the geometrical efficiency from 100% results from the fact 
that not all muons pass the upper plate vertically. By requiring a coincidence 
between the upper and the lower panel, the measured number of muons will 
therefore be less than the number of muons passing through the upper panel. 
The deviation of the geometrical efficiency was determined by simulating two 
surfaces with size of the scintillator panels with a separation of 30.5 cm. The 
start positions for the muons were randomly chosen on the upper panel. The 
azimuthal direction was sampled randomly in an interval between and 27r. 
The simulations were carried out with a lateral distribution according to equa- 
tion 2. The depth parameter h was varied o ver a range of ±10%. With a mea- 
sured density of (2.3±0.2) g/cm^ (|LAB9.5I) an efficiency of (88.5±0.2)% was 



obtaii ied. The vertical flux can be calculated using the approximating formula 
from (|MIY73l ) 



^{h,9) = <!){h,0)cos^-^'^e, (6) 
for a depth of 1526 m.w.e. the vertical flux is ^vert = (0.65 ± 0.04) x ^tot- 

3. 1 Backgrounds 



In the following section other sources of background than discussed in previous 
sections are estimated. In order to generate a coincidence pulse with an energy 
deposition of more than the threshold of 2.0 MeV one has to look at particles 
with higher energy. 



3.1.1 Neutrons from Muons 

Neutrons generated from muons traversing through the rock and salt are one 
candidate for the expected background in the WIPP. In his paper Bezrukov 



f)BEZ73^ estimates the amount of neutrons per muon generated at a depth of 



1500 meter water equivalent (m.w.e.) in hg/cm^ to 



iV„(^) = 7 10-^^. (7) 



With equation the empirical equation from 
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where $ has the units of (cm s sr and h represents the depth in hg/crn^, 
this calculates to a neutron flux of 3.36 10"'' kg~^s~^. With an assumed average 
cross section of 1 barn the attenuation length in salt is 21 cm. This calculates 
to a neutron flux of 1.6 10~^ cm~^s~^, or 8.6 10~^ cm~^ during the entire 
length of experiment of 6.17 days. Several Monte-Carlo simulations have been 
carried out to estimate the detection efficiency for neutrons in the existing 
detector geometry. The neutrons were randomly started from a sphere of 2 m 
radius surrounding the two panels into the whole solid angle. Several simula- 
tions with neutron energies ranging from 500 keV to 1 TeV were performed 
(see Fig. 8). According to the experimental settings, only such events were 
counted, where both panels had an energy deposition between 2 and 20 MeV. 
The most important interaction mechanism is inelastic scattering on hydro- 
gen and carbon in the scintillator material. Therefore the number of events 
depositing at least 2 MeV in each of the panels is reduced by at least 3 orders 
of magnitude, if the incident neutron energy falls below 4 MeV. Only neutron 
captures and other nuclear channels can than produce the required energy. 
The number of simulated events in Fig. 8 corresponds to 20 neutrons/cm^. 
Scaling this number to the expected number of neutrons during the experi- 
ment and assuming the highest efficiency, gives an upper limit for the number 
of neutron induced events of Nneutron = 2.6. Most of the neutrons will be 
moderated by the time they reach the detector. The assumption made above 
is, therefore, very conservative. 



3.1.2 Neutrons from U and Th 

Natural radioactivity is present in every geological layer of the earth. An 
abundance of uranium and thorium in the salt is able to generate neutrons 
via {a, n)-reactions and throug h spontan eous fission. The U and Th contents 
at WIPP have bee n mea sured ( WEB98[ l and are displayed in Table 3. With 



the numbers from (|FLQ88) for neutron production via (a, n)-reactions of U 



and Th in salt, who lists U(q,^„)=5.07 10"^ s"^g~^ for an uranium abundance of 
0.3 ppm and Th(Q „)=1.51 10"^ s^^g"^ for a thorium abundance of 2.06 ppm in 
salt, one can estimate the neutron production rate to U(q,„)=8.11 10~^ s~^g~^ 
for an uranium abundance of 0.048 ppm and Th(Q, ,„)=1.83 10"^ s^^g"^ for a 
thorium abundance of 0.25 ppm in salt (according to column "Used Values" 
in Table 3). Assuming again an attenuation length of 21 cm leads to an upper 
neutron flux of 1.28 10"^ s^^g"^. Since the highest alpha generated in the 
U/Th decay chain has ca. 8 MeV and the Q- values of the (n,a) reactions 
at Na or CI is at most -3 MeV, a maximum neutron energy of 5 MeV can 
be assumed. With the simulation results from Fig. 8 this leads to a neutron 
number in the panels of Nutk = 54 during the measurement. 
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3.1.3 Gamma Background 

There are two main sources of gamma background, first the decay of U/Th, 
and second the decay of ^°K. The highest level of gamma background comes 
from the decay of ^°K, which is abundant in the surrounding salt. The abun- 
dance of potassium in salt is 784 /ig/g (see Table 3) and corresponds to 
1.5 10^^ ^°K/g. A 1.46 MeV 7-ray will be emitted during the decay, with 
a probability of 10.7%. With an half-life of ti/2 = 1.28 10^ y this results in a 
decay rate of appro ximately 2.7 10~^ s~^g~^. With the macroscopic scattering 



cross section from (jBEROlh of cr=5.02 10^^ b in NaCl and the assumption of 
an exponential attenuation law, an attenuation length of 8.4 cm is computed. 
The attenuation length is defined as the length after which the intensity drops 
by a factor e. With a salt density of 2.3 g/civ? the amount of salt contribut- 
ing 7-rays can then be estimated to 19.3 g/cm^, which results in a 7-flux of 
$^ = 5.22 10^^ s^^cm"^. Since each gamma can deposit at most 1.46 MeV, 
a triple coincidence during the coincidence window of 100 ns in the detector 
is required to generate a signal inside the muon-cut region. According to the 
GEANT simulations, a 7-ray between 1 and 2 MeV emitted from a sphere 
of 2 m radius around the panels deposits energy with a probability less than 
1% in one of the detectors, resulting in a rate of 262 s~^ energy-depositing 
events. The resulting rate of triple coincidences with r = 100 ns and 

u = 262 s~^) is 1.80 10~^ s~^, or 0.1 events during the experiment, hence 
negligible. 



The number of gammas originating from (n, a) and (n, 7) reactions can be 
estimated using the above calculated neutron flux. Assuming that only one 7- 
ray above 5 MeV is produced per neutron, the flux is 1.28 10~^ s~^g~^. Gamma 
rays with energies up to 10 MeV have been simulated in order to estimate the 
efficiency of the muon panels including the coincidence requirement of at least 
2 MeV energy deposition per panel. The geometry and number of particles 
were the same as described above for the neutrons (see Table 4). Assuming 
that the detection efficiency was on average below the simulated efficiency for 
a 10 MeV 7-ray, at most of 32 counts are expected during the entire run of 
the experiment. 



3.2 The Muon-Flux 



The raw data have (5202±72) counts in the region of interest. After the correc- 
tion for electronic noise inside the region of interest (40±4) events and muon- 
induced events outside this region (62±6) events one obtains (5224±72 ± 7) 
counts. A conservative 10% uncertainty was assumed for each of the correc- 
tions, which were then quadratically combined (see Table 5). Only upper limits 
for the contribution from neutrons (55) and 7-rays (32) have been performed. 
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The uncertainty due to the detection efficiency is (±55) events. Since all the 
uncertainties are independent from each other, a quadratical addition of the 
different components can be performed. The final number of muon induced 
events in the muon panels is therefore (5224l;^Qg). As discussed above, the ge- 
ometrical efficiency, resulting from the fact that not all muons are vertically, 
has to be taken into account. The number of muons passing through the de- 
tector area is (5902lj23) events. The uncertainties of both, the surface area of 
the detector and the running time is well below the 1% level and can therefore 
be neglected. 

Now the muon flux can be calculated: 



" 532800 s X 23225 cm2 

= (4.77;°:?^) X 10-^ 1^ (9) 



This flux can now be converted to a vertical flux of 



Hz 

= (3.101°-°^) X 10-^ (10) 
cm^ sr 



converting to a meter water equivalence of (1585^6^) m. Me asureme nt s made in 
the p ast in dif ferent underground laboratories in the world ( CR,087 ). ( AMB95I ) 



and (|AND87[ ) show similar results for this depth. The shallow depth agrees 
within its error bar with a ffi through the previous experiments. Fig. 9 shows 
WIPP in comparison to other underground laboratories. With its average 
density of (2.3±0.2) g/cm^ and a depth of 655 m our measurement ffis very 
well with the predicted data. 



4 Conclusion 



The work described here has determined the vertical muon flux at the WIPP 
mine to ^vert = (3.10]'^o;o7) 10""^ s~^cm~2sr~^. This result agrees well with 
measurements at other underground laboratories and can now serve as an 
experimental basis for low-background experiments. 
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Table 1 

Cuts applied to determine event numbers in muon peak. 



Cut 


Cut 


Cut 


Events 


Count rate 


Number 


Upper 


Lower 






[Hz] 




Panel 


Panel 











E<2 MeV 


E<2 MeV 


5207± 72 


9.8±0.14 


xlO-3 


1 


— 


— 


11226± 106 


2.11±0.02 


xl0"2 


2 


E<2 MeV 




5881±77 


1.10±0.014 


xl0"2 


3 




E<2 MeV 


5350±73 


1.00±0.014 


xl0"2 


4 


E>2 MeV 




5345±73 


9.7±0.14 


xl0"3 


5 




E>2 MeV 


5876±77 


9.8±0.14 


xl0"3 


6 


E>2 MeV 


E>2 MeV 


5202±72 


9.9±0.14 


xlO-3 



Table 2 

Tail fits of the background and muon signal. 



Tail Name 


Fit 


Interval 


counts 


Monte- Carlo cut off 


0- 


-2.0 MeV 


62 


background upper panel 


1.2- 


-2.1 MeV 


20 


background lower panel 


1.0- 


-2.4 MeV 


40 



Table 3 

Natural radioactivity at the WIPP-underground. The abundance o f the ele ments 
was measured with mass spectrometry and gamma spectrometry by JwEBii). The 
Averages for the measurements are shown in the first two columns. The column 
labeled Used Value displays the value recommended by ( WEB98). The values for 
the columns in soil were typical measurement of regular ground dirt gathered in 
southern New Mexico. The las column displays the ratio of the typical soil value 
over the corresponding WIPP number. 





at WIPP 


Range in Soil 


Ratio 


Element 


Mass 

Spec. 

m 

9 


Gamma 
Spec. 

9 


Used 

Value 

M 

a 


low 

M 

a 


high 

M 

a 


typical 

M 

a 


Soil 
vs. 
WIPP 


Uranium 


0.048 


<0.37 


0.048 


0.5 


2.5 


1.5 


30 


Thorium 


0.08 


0.25 


0.25 


1.2 


3.7 


2.4 


10 


Potassium 


784 


182 


480 


500 


900 


700 


1.5 
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Table 4 

Number of events fulfilling the requirement of at least 2 MeV energy deposition per 
panel. 10^ gammas and neutrons were started from a sphere of 2 m radius around 
the muon panels. 



Enerffv 


Events 


5 MeV 





1 MeV 





2 MeV 





3 MeV 





4 MeV 


1 


5 MeV 


306 


6 MeV 


621 


7 MeV 


814 


8 MeV 


916 


10 MeV 


947 



Table 5 

The contributing parameter for the muon fiux. 

Nhvuc- of iiaramctor 



C<)u( ril)uli()n 



Runtime of the experiment 
Events in muon peak 
Contribution background to signal 
Muon induced events below threshold 
Neutrons from U and Th 
Neutrons from muons in rock 
Gamma rays 
Detector Efficiency 
Detector Efficiency (geometric) 
Conversion factor to vertical flux 



532800 seconds 
5224±72 Events 
-40±4 Events 
+62±6 Events 
-54 Events 
-2 Events 
-32 Events 
100 ± 1% 
84.25% ± 0.4% 
0.65 ± 0.04 
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Fig. 1. Setup of the scintillator panels. They were placed on top of each other with 
a distance of 30.5 cm. 



Panel 1 




Panel 2 



PMT 



777 



PMT 



777 



S. A. 



Disc 



coincidence 
unit 



S.A. 



Gate Delay 
Generator 



Fig. 2. Setup of the scintillator-clcctronics. The electronics was setup in a quad 
coincidence mode. The signal coming from the photo-multiplier-tube (PMT) was 
amplified by a Phillips Scientific fast amplifier (777) with a gain of 20. This output 
signal was split. One signal was used as input for a discriminator attached to a 
coincidence unit to form the strobe for the LeCroy AD811 ADC. The other output 
was run through a shaping amplifier (S. A.) into the input channels of the ADC. 



12 



"L 
"L 
"L 



"L 

I window 



Generator 



output 



"Too ^00 ^00 "400 ^00 ^00" 

Fig. 3. Timing of the scintillator Coincidence. 
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Fig. 4. Differential muon intensity versus muon momentum. The spectrum is an 
empirical fit of the momentum distribution of vertical muons at sea level. The 
spectrum has a flat maximum at about 0.5 GeV/c. 
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Upper Panel 




IS IS 20 
Eneigy [MeV] 



Lower Panel 




IS IS 20 
Eneigy [MeV] 



Fig. 5. Spectrum for the raw energy in the two scintihator panels. The energy 
plotted into a histogram in units of MeV, the Monte-Carlo simulation is shown 
dotted line. 
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15 20 25 

Energy lower panel [MeV] 



Fig. 6. Two dimensional spectrum for the cuts applied to the muons. The spectrum 
shows on the x-axis the deposited energy in MeV for the lower panel, on the y-axis 
the deposited energy for the upper panel. The black line displays the cut chosen 
for the muon number. The events marked by a cross represent the muon peak, the 
events marked by a dot are dominated by the random coincidence of the PMTs. 
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Energy lower panel [MeV] 

Fig. 7. Two dimensional spectrum for the cuts applied to the muon Monte-Carlo 
simulation. The spectrum shows on the x-axis the deposited energy in MeV for 
the lower panel, on the y-axis the deposited energy for the upper panel. The black 
line displays the cut chosen for the muon number. The events marked by a cross 
represent the muon peak, the events marked by a dot represent muons outside the 
cut region. The ratio of these muons in the simulation is 1.24 ± 0.01%. 
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Fig. 8. Neutron efficiency of tfie two panels as a function of primary neutron energy. 
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Fig. 9. Vertical muon flux for different underground experiments. Tfie errors are 
smaller than the data marker. 



18 



